Berry Phase Effects on Dynamics of Quasiparticles in a Superfluid with a Vortex 
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We study quasiparticle dynamics in a Bose-Einstein condensate with a vortex by following the 
center of mass motion of a Bogoliubov wavepacket, and find important Berry phase effects due to 
the background flow. We show that Berry phase invalidates the usual canonical relation between the 
mechanical momentum and position variables, leading to important modifications of quasiparticle 
statistics and thermodynamic properties of the condensates. Applying these results to a vortex 
in an infinite uniform superfluid, we find that the total transverse force acting on the vortex is 
proportional to the superfluid density. We propose an experimental setup to directly observe Berry 
phase effects through measuring local thermal atoms momentum distribution around a vortex. 
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Quasiparticles in a superfluid play an essential role in 
its thermodynamic as well as dynamic properties Q . Ow- 
ing to the realization of Bose-Einstein condensates (BEC) 
with ultra-cold atoms, much detailed investigation of the 
quasiparticles becomes possible |2j, |3j . In the past several 
years, extensive experimental and theoretical work has 
been devoted to the study of quasiparticle dynamics in 
static BECs |U d S- E3| - while quasiparticles in a flowing 
condensate have not received much attention 0. 

In this Letter, we study the quasiparticle dynamics in 
a BEC superfluid with a vortex, a typical flowing con- 
densate. We adopt a semiclassical approach by following 
the center of mass motion of a Bogoliubov wavepacket, 
finding important Berry phase effects Q due to the 
background flow. Unlike the case with a static BEC, 
the quasiparticle position and momentum are no longer 
canonical variables , and there are also important modifi- 
cations to the quasiparticle statistics and thermodynamic 
properties of the condensates. We use these results to cal- 
culate the momentum circulation around a vortex in an 
infinite superfluid and find the total transverse force act- 
ing on the vortex is proportional to the superfluid density 
which rules out the existence of the Iordanskii force 
that is supposed to originate in the asymmetric scattering 
of quasiparticles by the vortex [ToL fill fl2T ] . Finally, we 
propose an experimental setup to directly observe Berry 
phase effects through measuring local thermal atoms mo- 
mentum distribution around a vortex. 

Wavepacket dynamics — For simplicity, we present the 
theory here for a single vortex, although it is applicable 
to a general flowing condensate such as vortex lattices. 
The dynamics of the condensate are described by the 
Gross-Pitaveskii equation [LH ] 
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atom coupling constant, and V (r) is the trapping po- 
tential. We have also included an angular momentum 
term to take into account the possibility of using a ro- 
tating reference frame, with f2 being its angular fre- 
quency. The condensate wavefunction may be written 
as ip = \J n (r)e la ^e~ ltJ/t , where fj, is the chemical po- 
tential, n (r) and a (r) are the vortex density and phase 
respectively. The gradient of phase yields local velocity 
of the vortex v s = Vet (r) . 

Small deviation 8ipe~ ll - lt (elementary excitations) from 
the condensate wavefunction tpo is described by the Bo- 
goliubov equation 
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The wavefunction $ has two components, u and v, which 
are related to 5vj through Sijj = ue~ %ut + v*e iuJt , where 
u! is the quasiparticle energy. The entries of the matrix 
operator are given by H± = p 2 /2 + 2gn(r) + V (r) — 
/i =p (r x p), and H2 (r) = gn (r). Because of the Pauli 
matrix a z on the left hand side, the Bogoliubov equation 
is nonhermitian. 

We now consider a quasiparticle wavepacket centered 
at r c with its spread small compared to the length scale 
of the slowly varying potentials (including trap poten- 
tial V (r), condensate wavefunction tpQ and terms re- 
lated to n). The dynamics of the wavepacket is ap- 
proximately governed by the local Bogoliubov operator 
Q c = Q(p, r = r c ) plus its linear gradient correction. 
The local Bogoliubov operator has plane wave eigenstates 
Q c e lqr |0(q,r c )) = uj c a z e % * T \<j> (q, r c )}, where q is the 
wavevector. The amplitude satisfies a simple 2x2 matrix 
equation, which can be solved easily to yield 



(q,r c ) = (ff 2 (q)-ff 2 2 
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where for simpler notation we have used units such that 
h and the atomic mass m are both unity, g is the inter- 



for the local quasiparticle energy, where H\ (q) = 
q 2 /2 + 2. 9 n(r c ) + V (r c ) - and |0(q,r c )) = 
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tude of the local eigenstate, where C — i h 1 +h 2 
note that the wavevector q and the wavepacket center 
r c enter the local quasiparticle energy and wavefunction 
parametrically. The wavefunction is normalized in the 
sense that (</> (q, r c )| <r z \<j> (q, r c )) = 1. We have also cho- 
sen the phase of the wavefunction such that it is smooth 
and single valued in the parameters (q, r c ). We shall see 
that the parametric dependence of the eigenstates on the 
center position of the wavepacket will manifest as Berry- 
phase terms in the equations of motion. 

We now turn our attention to the wavepacket itself, 
which is to be constructed out of these eigenstates as 



|3>= J dV(q,t)e^ r |<Mq, 
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where the superposition amplitude a (q, t) may be taken 
as a Gaussian in q. The normalization is taken to be 
|<f>) = J d 3 q \a (q,i)| 2 = 1. We assume that the 

|2 



Gaussian is centered at q c = j d 3 q \ a (q, t)\ 2 q and has a 
width narrow compared to momentum scales of the en- 
ergy dispersion and of the eigenstates. Microscopic calcu- 
lation shows that q c = ($| a z p |$) = J Sip*pSipdr, imply- 
ing that it represents the mechanical momentum of the 
quasiparticle. To be self consistent, the wavepacket must 
yield the preassigned center position r c = ($|cr z r|$). 

The dynamics of the quasiparticle can be derived from 
a time-dependent variational principle for the Bogoliubov 
equation, where the action (defined as time integral of the 
Lagragian) is extremized with respect to the quasiparticle 
wavefunction. Here the Lagrangian is given by 
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and the wavepacket (0J is chosen as the variational wave- 
function with time-dependent parameters r c and q c . We 
use d/dt to mean the derivative with respect to the time 
dependence of the wavefunction explicitly or implicitly 
through r c and q c . Under the previously discussed con- 
ditions that the wavepacket is narrow both in position 
and momentum spaces, the Lagrangian JSJ can be eval- 
uated as a function of variational parameters r c and q c , 
and their time derivatives, independent of the width and 
shape of the wavepacket in position or momentum. 

Quasiparticle energy and Berry phase — The first term 
in Lagrangian (JSJ corresponds to the total energy of the 
quasiparticle wavepacket, which is found to be 

lu = [Hi (q c ) - H 2 ] 1/2 - q c • (fi x r c ) + (l - p 2 ) q c • v s . 

(6) 

where v s = Va (r c ) is the local velocity of the vortex. 
The first two terms are the local quasiparticle energy 
Eq.© with momentum q c , stemming from the expecta- 
tion value ($| Q c $) of the local Bogoliubov operator. 



The last term is the correction of the wavepacket en- 
ergy originating from the linear gradient expansion of 

the Bogoliubov operator AQ = \ {(x — r c ) ■ 0^ + c.cY 

and p = ($|$) is the total atomic mass contained in the 
quasiparticle wavepacket. 

To understand the quasiparticle energy expression, it 
will be instructive to consider the simpler situation of 
a uniformly flowing superfiuid of velocity v s . If po is 
the quasiparticle momentum in the reference frame where 
the superfiuid is static, then we have the relationship 
qc = Po + pv s . In terms of po, the quasiparticle energy 
© can be expressed (in the limit of small v s ) as w = 

1/2 

e (Po) + Po ■ v s , where e (p ) = [Hi (p ) -iff] = 

1 /2 

[p'o (.Po/4 + gn)] is the energy dispersion in a static 
superfiuid. This is just the standard Landau formula 
obtained using Galilean transformation 

The evaluation of the second term of the Lagrangian 
(0 is similar to that in Ref. 0| > which yields 
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Here the term q c ■ f c comes from the time-dependence 
of the superposition amplitude a (q, t) in the wavepacket 
Q, while the vector potential A = i (<j>\ a z \dcf>/dr c ) = 
— (p— l)Va (r c ) is obtained from the position dependence 
of the two-component local eigenvector, which is similar 
to the case of a spin in a position-dependent Zecman field 
[8j . The vector potential, also called Berry connection, 
arises for a vortex because of the non-zero local veloc- 
ity. Its line integral over a path C gives a Berry phase 
T (C) = § c dr c ■ A = — § c (p — 1) da (r c ) of the eigenvec- 
tor []J| . For instance, the accumulated Berry phase for a 
quasiparticle moving around a vortex in an infinite uni- 
form superfiuid is — 2ir (p — 1), where p represents the to- 
tal atom number (notice that the atom mass m has been 
taken as unity) in the quasiparticle wavepacket. Apart 
from a non-essential 2tt, it means that the phase acquired 
by moving each atom around a vortex is 2tt, which agrees 
with Feynman's argument |l6j . 

We note that the vector potential vanishes at very high 
momenta, where p — > 1. Likewise, the energy correction 
due to superfiuid flow (the last term of Eq.©) vanishes 
in this limit. This is quite reasonable, because the quasi- 
particle becomes a free particle and is decoupled from the 
condensate at high momenta. The condensate has influ- 
ence on the quasiparticle only at low and intermediate 
momenta. 

Semiclassical dynamics — Following the standard pro- 
cedure of analytical mechanics, we may introduce the 
canonical momentum k c = dL/ di c — q c + A conjugate 
to the coordinate vector r c . We see that the vector poten- 
tial makes the canonical and mechanical momenta differ- 
ent. The energy lu = k c - f c — L becomes the Hamiltonian 
and the equations of motion are 

• >L: <h: (8) 
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where the quasiparticle energy ui is now recasted as 

a function of k c , lu = [iff (k c ) — flf] 2 — (k c — A) • 
(f2 x r c ) + k c • A. We see that the vector potential A 
modifies the energy expression in this canonical formula- 
tion. 

In terms of the physical quantity of mechanical mo- 
mentum q c , the quasiparticle equations of motion are 
dramatically altered: 



duj dp 
dq c dq c 
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duj . f dp \ ( . dp 

q c = -„ r c x — x v s + q c • — | v.,. 

or c \dr c ) \ dq c 

They follow directly from (JHJ by the transformation 
k c = q c + (1 — p) v s . Here oj is the original energy ex- 
pression ® . We see that additional terms depending on 
superfluid velocity and total atom mass in the quasipar- 
ticle wavepacket appear in both equations, therefore the 
equations of motion are no longer of canonical form. This 
a result very different from the case of static superfluid 
studied in 

In statistical treatment of quasiparticles, an important 
semiclassical quantity is the density of states of phase 
space. For the canonical variables k c and r c , this should 
be taken as the constant 1/ (2ir) 3 . Because of the vector 
potential, the position r c and mechanical momentum q c 
are no longer canonical variables, which leads to a mod- 
ification of the density of states |l7| 
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where I is the 3x3 unit matrix. 

The thermodynamics properties of condensates are 
determined by the physical quantities of quasiparticles 
such as energy and density of states, therefore they are 
strongly affected by Berry phase. For instance, the ther- 
mal depletion of the condensate density is given by 



Sn (r c ,T) = J d 3 q c D(r c ,q c )p/ (t 
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according to our theory, which is different from that for a 
static condensate because of the modifications of energy 
u> and density of state D (r c , q c ). Notice that the inte- 
gration of this distribution over position space yields the 
total number of thermal atoms. 

Transverse force on a vortex — The general form of 
the transverse force per unit length acting on a moving 
vortex with velocity in an infinite uniform superfluid 
can be written as F = k (C) z x v l , where k (C) is the 
momentum circulation along a path far away from the 
core of the vortex, and we assume the normal fluid ve- 
locity is zero |l2|. In the past several decades, there 
has been some controversy about the expression of k (C), 
which was argued to be either 2Trhn tot /m or 2Trhn s /m 



through different approaches 0,0,^3- Here 2-Kn/m is 
the quantum of circulation, n to t is the total mass density, 
and n s is the superfluid density . 

According to our theory, the total momen- 
tum circulation k (C) contains two contributions: 
2irh (ntot ~ (r c , T)) /m from the condensate and 
K\ (C) = § c W • dr c from quasiparticles. W is the 
momentum density of quasiparticles under equilibrium 
Bose-Einstein distribution and is obtained by summing 
over the mechanical momenta from all states, 
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where T is the temperature, fcg is the Boltzman constant, 
D(r c ,q c ) = ^(l-^frcxq,),), and 1/ (2tt) 3 

are the density of states for the non-canonical momen- 
tum q c and canonical momentum k c respectively. 

Since we choose the circular contour C to be far away 
from the core of the vortex, we can expand the Bose- 
Einstein distribution / (k c ) w fo (k c )— A/ with /o (k c ) = 
-. ,} T : , and A/ = . e °' ^ — (i-p)k^-ve c c ons jj_ 

e »e/»B' - 1' •' ( e ^ c /k B T_ 1 y k B l 

ering a closed orbit with fixed radius R and taking the 
limit of R — > oo, we find the thermal depletion of the con- 
densate Sn(r c ,T) — ^2 Jo" (k c ) pk^dk c , and quasi- 
particle momentum circulation 



Kl (C) 
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The first term of K\ (C) stems from the difference be- 
tween the canonical momentum k c and mechanical mo- 
mentum q c , and cancels with — 2TrhSn (r c , T) /m. The 
second term of K\ (C) comes from the deviation A/ in 
the distribution function, and is found to be 2nhp n /m, 

9 2 (k T) 4 

where p n — -jg- ^ J is the normal fluid density. Sum- 
ming up the contributions from both condensate and 
quasiparticles, we find the total momentum circulation 
k(C) — 2irh(n to t — Pn) /m = 2itTin s /m, and the trans- 
verse force per unit length F = 2irft,n s z x v^/m. This 
form of transverse force agrees with that obtained by 
Thoulcss et. al. using general properties of superfluid 
order 0, and rules out the existence of the Iordanskii 
force that is proportional to the normal fluid density p n , 
and is supposed to originate in the asymmetric scattering 
of quasiparticles by the vortex [lfj, [ll], Q3 • 

Experimental observation - The above discussions 
show that the modifications of mechanical momentum, 
energy dispersion and density of states of quasiparticles 
due to Berry phase affect the thermodynamics of atoms. 
Our prediction of a transverse force proportional T can 
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Figure 1: Thermal-atom momentum distribution around a 
trapped vortex. Length scales are in units of I — 0.73 /im for 
a Rb 87 BEC in a trap specified in the text, (a) A schematic 
plot of experimental geometry. Measurement of thermal-atom 
momentum distribution is to be made in a region S at 100 
unit lengths away from the vortex center with a radius of 
20 unit lengths. This region contains about 9740 conden- 
sate and 1650 thermal atoms, (b) Condensate atom density 
(multiplied by the interaction strength g) in the presence of a 
trapped vortex, (c) Density of states for quasiparticles in the 
measurement region S as a function of quasiparticle momen- 
tum q c along the circulation flow, (d) Momentum distribution 
of thermal atoms in the measurement region S. Dotted and 
dashed lines correspond to our Eq. I14H and the naive ex- 
pression J s dr c p/ (exp ( " c ^r Ve ) — l): respectively. For ref- 
erence, the momentum of condensate atoms is indicated by 
the vertical solid line. 



in principle be measured in experiments. In the follow- 
ing, we propose an experiment for a direct observation of 
the local thermal atom momentum distribution around 
a vortex, in which the effect of Berry phase can be seen 
clearly. 

Fig. 1(a) is a schematic plot of experimental geometry. 

The length unit is / = A / — , h f , chosen for convenience. 

It corresponds to I « 0.73 /im for a BEC of TV = 5 x 10 5 
Rb 87 atoms in a quasi-two dimensional magnetic trap 
|l8| with axial and radial trapping frequencies oj z — 2ir x 
800 Hz, oj r = 2tt x 2 Hz, in which a single vortex is 
created at a rotation frequency of fl = 0Au> r [l9]. The 
condensate density profile is shown in Fig. 1(b), where one 
can see that the observation region S is chosen relatively 
far away from the vortex core, which was to make sure 
that the semiclassical approximation remains valid. We 
will consider a temperature of T = 2fi/kB ~ 21.1 nK 
which is about a third of the BEC transition temperature 
T c k 52.8 nK for this system. 

The atoms are supposed to be in a hyperfme state (|1)). 
By focusing two copropagating Raman beams at region 
S one may drive both the thermal and condensate atoms 



there to another hyperfme state (|2)). The remaining 
atoms in state 1 1) are to be rapidly expelled from the trap 
by applying radio-frequency radiation that flips them to 
an anti-trapped hyperfme state. Finally, the time of flight 
measurement 20] of atoms in state |2) can then deter- 
mine their momentum distribution, which, according to 
our theory, is given by 

Sn (p c , T) = f d 2 r c D (r c , q c ) p/ (f' k ° T - l) . (14) 

In Fig. 1(c), we see that there are more quasiparticles 
moving along the circulating superfluid flow than in 
the opposite direction and the deviation from the naive 
value of the density of states 1/ (27r) 2 is substantial. In 
Fig. 1 (d) , we find that the momentum distribution is quite 
different from the commonly believed one that does not 
take into count the Berry phase effects. 
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